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ABSTRACT 

We construct a family of simple analytical models of galaxy clusters at the present 
epoch and compare its predictions with observational data. We explore two processes 
that break the self-similarity of galaxy clusters: systematic variation in the dark mat- 
ter halo concentration and energy injection into the intracluster gas, through their 
effects on the observed properties of galaxy clusters. Three observed relations between 
cluster properties and temperature are employed to constrain the model; mass, slope 
of gas density profile (/?) and luminosity. The slope of the mass-temperature relation 
is found to be reproduced by our model when the observed variation in concentration 
is included, raising the logarithmic slope from the self-similar prediction of 1.5, to 
that of the observed relation, ~ 2. Heating of the intracluster gas is observed to have 
little effect on the mass-temperature relation. The mean trend in the /^-temperature 
relation is reproduced by energy injection in the range 0.5-0.75 keV per particle, while 
concentration variation is found to have only a small effect on this relation. Excess en- 
ergies calculated for individual systems from the /3-temperature relation suggest that 
the lowest mass systems may have excess energies that are biased to lower values by 
selection effects. The observed properties of the luminosity-temperature relation are 
reproduced by the combined effects of excess energy and a trend in the dark matter 
concentration. At high masses the observed variation in dark matter concentration 
results a logarithmic slope of ~ 2.7 compared to recent observations in the range 2.6- 
2.9, whilst the observed steepening of the relation in galaxy groups is predicted by 
the model when heating in the range 0.5-0.75 keV per particle is included. Hence a 
combination of energy injection and dark matter concentration variation appears able 
to account for the mean trends in all the observed relations. Scatter in the energy 
injection and concentration may account for a large proportion of the scatter in the 
observed relations. 
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1 INTRODUCTION 



In the ' liiei ai ihieal clustering model for — the foi illation of 

structure in the universe, small scale perturbations collapse 
into virialized objects and then cluster together to form suc- 
cessively larger virialized structures. Simple physical models 
of gravitational collapse along with shock heating of the gas 
would suggest that these virialized structure s should be ap- 



prox imately scaled versions of one another (Navarro et al 
1995). In the absence of other physical process such as heat- 



ing or c ooiing, the gas and dark matter halos will be almost 
self-similar, it is possible tor energy to be transferred be- 
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tween the dark matter and the gas during mergers (Pearce 



et al. 1994 ) which is believed to be responsible for the cores 



observed in cluster gas density profiles. However this is not 
usually assumed to act in a way which would break self- 
similarity. 

This expected self-similarity has already been observed 
to be broken in several respects. The surface brightness 
profiles of galaxy clusters have been observed to flatten in 



low m ass systems (Ponman et al. 1999; Helsdon & Ponman 



2000b) and their gas density, entropy and energy have been 
shown not to be self-similar ( |Llo y d-Davies et al. 200c| ). The 
luminosity-temperature relation for galaxy clusters is also 
observed deviate from its expected self- s imilar behaviour 
(|Edgc fc Stewart 199 1|; |Markevitch 199^; |ficlsdon fc Pon ' 
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man 2000a). These effects are usually attributed to energy 
injecti on into the intracluster medium (ICM) from galaxy 
winds (Lloyd-Davies et al. 2000). However there is some con- 
troversy over whether galax y winds can provide the energy 
nece ssary to heat the IC M (Valageas & Silk 1999; Wu et al 



200C; Eower et al. 2001). An alternative often promulgated 
is that AGN are responsible for the injected energy. How- 



ever wqile the amount of energy available from AGN is ver 



rv 



which results in flattening of the gas density profile, and 
systematic variation in the dark matter halo concentration. 
In order to put constraints on these processes one useful 
approach is to construct models of galaxy clusters and then 
compare the predictions of the models with observed proper- 
ties of clusters. A similar a pproach has previously been used 
by a numb er of authors JCavalicrc et al. 1997^ Wu et al. 
1991, iBoU ICavaliere et al. 19991; palogh et al. 1999| ; [Tozzl 



the ICM. AGN have emitted large amounts of energy in the 
form of electromagnetic radiation, especially in the ultra- 
violet, but this will not heat gas to the required tempera- 
tures. Mechanical energy transfer through AGN jets seems 
the only possible option at present, but it is not clear how 
much energy is available in this form or whether it can be 
transferred efficiently into the ICM. It shou ld be noted that 
preliminary Chandra results ( ^abian 2001 ) do not suggest 
that widespread heating of the ICM by AGN is occurring at 
low redshift. 

The amount of energy injection inferred from the ob- 



large, ii is not entirely clear how they could transfer it into & Norman 2001; power*et*'al. 2001 ) to explore the effects of 



served breaking of similarity is also controversial. Wu et al 



(2000|)j|uggest a value of 1.8-3.0 keV per particle is needed to 
explain the steepening of the luminosity-temperature rela- 
tion from the self-similar prediction of ~ T 2 to the observed 



value of ~ T . However Lloyd-Davies et al. (200C) have mea- 



sured a mean excess energy of 0.44 keV per particle in a sam- 
ple of 20 galaxy clusters and groups using analytical models 
fitted to X-ray spectral images. There would also seem to 
be a theoretical argument against energy injection of much 
greater than 1 keV per particle since it is difficult to see 
how galaxy groups with mean temperatures less than 1 keV 
could have X-ray h aloes, whereas a number of such systems 
are known to exist ( Ponman et al. 19961; Pelsdon & Ponman 



2000b|] It therefore seems unlikely that energy injection into 
the the ICM can be the explanation for the steepness of the 
cluster luminosity-temperature relation or that any energy 
injection can greatly exceed 1 keV per particle. 

Breaking of simple self-similarity in the dark matter 
constituents of galaxy clusters has also been observed. The 
dark matter concentration of galaxy clusters systematically 



decreases with i ncreasing cluster mass (Wu fc Xue 2000 
|Lloyd-Davies fc Ponman 2002[ ). This effect 



Sato et al. 2000 



has been observed fo r some years in nu merical simulations 



( Navarro et al. 1997 ; Moore et al. 1999| ) . It is generally in- 



terpreted as due systematic variations in the formation time 
with system mass. This is predicted by the hierarchical clus- 
tering model of structure formation, with low mass systems 
collapsing first and then clustering together to form pro- 
gressively larger structures. The exact process by which this 
affects the dark matter concentration is a matter of some 
uncertainty, but an interesting hypothesis is that halos all 
form with the same concentration, and subsequent accre- 
tion increases the virial radius and therefore the concentra- 



tion (Salvador-Sole et al. 1998). Therefore halos that form 
at high redshift and therefore have had a long time to ac- 
crete material will have the highest concentrations. As the 
gravitational potential affects the temperature structure of 
the ICM, variation of the dark matter concentration would 
be expected to have an impact on the observed properties 
of the ICM. 

We therefore have two processes which are observed to 
break self-similarity in galaxy clusters. Heating of the ICM 



energy injection. However these models have in general been 
driven from a theoretical perspective and have also tended 
to concentrate their efforts on predicting the evolution of 
cluster properties with redshift. Our aim is to construct a 
model of galaxy clusters, base on empirical results wherever 
possible, which will predict their properties at zero redshift. 
We will then use this model to investigate in detail how pro- 
cesses that break self-similarity, such as injection of energy 
into the ICM, affect the systems' observed properties. Com- 
parison of the model predictions with observed data should 
then allow constraints to be placed on the various process 
involved in breaking self-similarity. 



2 CLUSTER MODEL 

In order to construct useful models of galaxy clusters to 
be compared with observations, a number of simplifying as- 
sumptions need to be made. In this case we will use only one- 
dimensional spherically symmetric models containing only 
dark matter and hot gas in hydrostatic equilibrium. While 
these approximations are crude, it should be noted that most 
analyses of observations of galaxy clusters involve the as- 
sumption of spherical symmetry. The approximation is fairly 
good for relaxed systems but it should be born in mind when 
considering the results that it is not a very good approxi- 
mation to morphologically disturbed systems. In our model 
we make no distinction between dark matter and galaxies. 
However to a first approximation galaxies can be considered 
to be collisionless and while dynamical friction may be have 
some effect on the galaxies in low mass systems, it is unlikely 
that this will significantly affect the gravitational potential 
given that the galaxies make up only a small fraction of 
the mass of the system. Our approach is not to model the 
evolution of galaxy clusters to the state we observe them 
at present, but to model a variety of possible end points of 
cluster evolution and then compare them against observed 
clusters. 

The primary parameter in the cluster model is the clus- 
ter mass. There are a number of possible ways the cluster 
mass could be defined. In general the mass within the region 
around the cluster centre that is in virial equilibrium is con- 
sidered. The overdensity of this region with respect to the 
critical density, p cr it, can be calculated for the collap se of a 
spherical top-hat density perturbation (Peebles 1980). For a 
critical density universe this overdensity is 187r (~ 178) but 
can be signific antly smaller for other cosmologies (Bryan & 
Norman 1998). An overdensity value of 200 is often used to 



define the outer boundary of a system as this is smaller than 
the virial radius for all reasonable cosmologies. In order for 
our results to be easily comparable with observations this 
overdensity was used. This results in a relationship between 
the mass, M200, and the radius, -R200, of 
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-R200 = 



3M 20 o 

800ivp cri 



(1) 



to be anti-correlated with cluster mass in numerical simula- 



tions (Navarro et al. 1997) and this has been o bserved using 



X-ray data (Wu fc Xue 2000; 



Sato et al. 200C 



Lloyd-Daviet 



which depends only on the critical density, p cr j t . Throughout fc Ponman 2002j). It is possible for us to take account of this 



the modeling we adopt Uo=t>U km s Mpc and go = 0.5. 



2.1 ICM density profile 

In our model the gas density of the ICM is represented by 
the usual parameterization of the form, 



variation in concentration using an analytical relation, for 
instance the relation 



34.9 



( Macro 
^1O 13 M 



(4) 



is measured by Lloyd-Davies fc Ponman (2002| ) for a sam- 



p[r) = p(0) 



1 + 



(0 



(2) 



where r c is the core radius, (3 is the density index an d p(0) 
is the density normalization. Jones fc Forman (1984 ) have 
found this to be a good representation of the structure of the 
ICM in galaxy clusters. Vikhlinin et al. (199£ ), using the best 
quality ROSAT data, agree with this, although they find 
marginal evidence for steepening of the surface brightness 
profile at large radius. The fiducial value for the j3 parameter 
in our model was the canonical value of §. A fiducial gas 



pie of 20 galaxy clusters and groups. It should be noted 
that this relation is a rough characterization of a trend that 
shows considerable scatter. It is also possible to use the re- 
sults of numerical simulations to define the dark halo con- 



centration. The simulations of Navarro et al. (1997) predict 
concentration variation for a ACDM cosmology that can be 
approximated by the relation 



c = 13.2 



/ Macro 



(5) 



fraction within i?aoo of 0.2, as observed by Lloyd-Davies & 



density profile 



Ponman (2002), was used to set the normalisation of the gas 



over the mass range of galaxy clusters and groups. This 
predicted relation is consid e rably flatter than the those ob- 
served by |Wu fc Xue (2000| ); |Sato et al. (200C| ); [Lloyd-Davies 



& Ponman (200i ). We will investigate the use of both rela- 



To define the core radius of the gas densi ty profile, r c , 
we use the value of 7 percent of i?aoo found by Lloyd-Daviet 



fc Ponman (2002) for systems where a cooling flow does not 
obscure the core. Selecting a value for the normalization of 
the gas density profile, p(0), is not a trivial problem. One 
possibility is to normalize the profile such that the gas mass 
within i?200 is f gas M2oo where f gas is the gas fraction of 
the cluster and a parameter to be specified. However this 
method may be unphysical as it does not allow gas to be 
pushed outside 7?aoo when the specific energy of the gas 
is raised by a large amount. The normalisation of the gas 
density is therefore a boundary condition for the problem 
which needs careful consideration. 



2.2 Dark matter density profile and gravitational 
potential 

In our model no distinction between dark matter and stars 
is made, since both are expected to be collisionless on a large 
scale and neither will contribute significant emission to the 
X-ray data we will be comparing the model against. The 
dark matter distribution in our model is re presented by a 



tions. The normalization of the of the dark matter profile, 
Pdm, is set so that the cluster contains a dark matter mass 
OI (1 _ fgas)M 2 oo within i? 2 oo- 

With the gas and dark matter profiles defined it is then 
possible to calculate the cluster temperature profile assum- 
ing hydrostatic equilibrium. For hydrostatic equilibrium and 
spherical symmetry, the equation 



M(r) = 



T(r)r 



dlnp dlnT 



dlnr dlnr 



(6) 



profi le derived from numerical simulations (Navarro et al 



1995) of the form, 



PDM(r) = Pdm 



[x(1 + i 



(3) 



where x — r/r s and r s is a scale radius . It should be noted 
that it has recently been suggested (Moore et al. 199!:) 
that higher resolution simulations produce dark matter pro- 
files with steeper central cusps, with asymptotic slopes of 



about 1.5 rather than 1 in the profile of (Navarro et al 



1995). However this effect is only important at small radii 



and as m any observational res ults are derived using the 
profile of ( Navarro et al. 1995 ) it is better suited to our 
purposes. It is possible to parameterize the concentration 
of the dark matter profile by th e concentration parameter, 
c = raoo/r s ( Navarro et al. 1997). This concentration is seen 



is satisfied (Fabricant et al. 1984). Given the mass and gas 
density profiles of a cluster the temperature profile is defined 
if the temperature is specified at some point. This introduces 
a second boundary condition to the problem that we must 
consider. 



2.3 Effects of ICM cooling 

In our model we do not explicitly model the effects of ra- 
diative cooling on the ICM for several good reasons. The 
amount of cooling in a particular cluster is will be related to 
the amount of time it has been left undisturbed by a major 
merger. This is an essentially random factor which varies 
from system to system, depending on its particular history. 
In contrast, our model aims to represent the present state 
of clusters, and to avoid the complications and uncertainties 
associated with cluster evolution. Secondly while the physics 
associated with radiative losses from ion-electron interac- 
tions in the ICM is well understood, the macroscopic effects 
of this cooling are not. The cooling is thought to result in 
a highly multiphase structure in the ICM but insufficient 
details are known to accurately model it. Thirdly, current 
observations do not provide enough information to provide 
a secure empirically based model of cooling in clusters. 

We therefore adopt an approach of modeling the clus- 
ters without any cooling included and then attempting to 
take account of this when comparisons are made with ob- 
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servational data. In many cases it is possible to use observa- 
tional data that has been corrected for the effects of cooling. 
This should allow a reasonably unbiased comparison of the 
model predictions with the data. In some cases, especially 
for galaxy groups cooling corrected data is not available. 
In these cases we will attempt to estimate what the likely 
effects of this are on our results using empirical data on 
correction for cooling. While this is far from ideal we be- 
lieve that attempting to model the effects of cooling would 
introduce at least a similar amount of uncertainty into our 
results, given our present state of knowledge of cooling flows. 



2.4 Boundary conditions 

For a cluster that has the specific energy of the gas raised 
above the default value by flattening its gas density profile, 
there are two boundary conditions that must be specified. 
A normalisation for the gas density profile and a normalisa- 
tion for the temperature profile. The default clusters need 
only one boundary condition, the temperature profile nor- 
malization, as the gas density normalisation is set by the 
gas fraction within the virial radius. However for the clus- 
ters with raise specific energy if the model is to allow the gas 
fraction to vary another constraint must be found. The most 
physically justifiable places to set the boundary conditions 
are at the centre of the systems and at the shock radius. 
However the position of the shock will be dependent on the 
amount of time that has elapsed since the system formed 
and also on the extent to which the infalling gas has been 
previously heated. Simulatio ns suggest that the shock radius 
should occur at 1 — 1.5-R„j r (Knight & Ponman 1997; Tozzi 



fc Norman 2001 ) and that preheati ng can result in the shoc k 
propagating out as far as 2.5R v i r Ibzzi & Norman (2001), 



since higher entropy gas has a higher sound speed. 

One possible way of setting t he density boundary con- 
dition is to use the observation of Lloyd-Davies fc Ponman 



(2002 ) that the gas density extrapolated to -R200 appears to 



converge to an approximately constant value for systems 
over a wide range of system masses. It should be noted 
that this result relies on extrapolating models fitted to data 
within a radius smaller than -R200 and is not a direct mea- 
surement of the density at i?200- The gas density at -R200 in 
the raised specific energy case can therefore be fixed at the 
density in the default case. This approach has the advantage 
of being extremely simple and has a least some justification 
from observations. We will therefore use this condition for 
the main results in the paper but wil l inv estigate the effect 



of alternative assumptions in Section 4.5. 

A constraint is also needed to set the temperature 
boundary condition for the model. Unfortunately observed 
temperature profiles of galaxy clusters are much less well 
constrained than gas density profiles and therefore there are 
not much in the way of observational constraints on the tem- 
perature at .R200- Therefore it is necessary to resort to theo- 
retical constraints on the temperature profile. A constraint 
Bower et al. (2001 ) derived from the numerical simulations 
of |Ekc ct al. (1998 ) and |Frenk et al. (200C| ) is that the tem- 
perature at the virial radius is approximately 0.5T„i r , where 
T V i r is defined by the equation, 



The mean mass per particle [i, is taken to be 0.6 amu. 



2.5 Energy calculation 



Observational studies (Lloyd-Davies & Ponman 2002; Hels- 



don & Ponman 2000b) suggest that the ICM density profiles 



in low mass systems are flattened compared to the profiles 
of high mass systems and this is observed as a reduction 
of the f3 parameter in low mass systems. There is little evi- 
dence to suggest, although it cannot be ruled out at present, 
that there is any significant effect on the core radius of the 
gas density profile. Figure |l| shows the core radii of the 8 
galaxy cluster and groups in the sample of Lloyd-Davies & 
Ponman (2002| ) with reliable core radii measurements, plot- 
ted against emission weighted temperature. The solid line 
shows the relation predicted by the model between core ra- 
dius and emission weighted temperature when the core ra- 
dius is a constant fraction, 7 percent, of i?200- It can be seen 
that while the constraints of the data are not very strong, 
a model where the core radius is a fixed fraction of -R200 
appears to be quite consistent with the data. Since observa- 
tional evidence strongly favo urs variation in the parameter 
rather than the core radius ( Lloyd-Davies fc Ponman 2002 



Helsdon fc Ponman 2000b ; Mohr et al. 1999), this approach 



therefore seems the most justifiable to take in modeling vari- 
ation in the gas specific energy. 

Our model therefore allows the gas specific energy to 
be raised by reducing the value of the (3 parameter and so 
flattening the gas density profile. In order to quantify differ- 
ence in gas specific energy from the default case, a prescrip- 
tion is needed for how the total energy of the ICM is to be 
measured. It is possible given the temperature and density 
structure of the ICM, and also the dark matter distribution, 
to calculate the thermal and gravitation al potential of the 



ICM. However it has been pointed out by Lloyd-Davies et al , 
(2000|) that care must be taken in selecting the regions over 



which to calculate the energy as if the gas distribution is 
changed the energy of the ICM within a fixed radius will 
not be comparing the same mass of gas. To calculate the 
difference in the energy of the ICM between raised specific 
energy and default models we therefore calculate the differ- 
ence between the total energy of the gas within -R200 in the 
raised specific energy model and the total energy of the same 
mass of gas in the default model, which will be contained 
within a smaller radius. Altering the value of j3 affects both 
the gravitational potential and thermal energy of the gas, 
since the gas temperature profile will be modified if the gas 
density profile is changed. Using this prescription it is pos- 
sible to alter the value of (3 until the desired difference in 
specific energy is ach ieved. This approa ch is similar to that 
used in the model of Bower et al. (2001). 



T v i r — 



2kR 20 



(7) 



3 OBSERVATIONAL DATA 

In order to constrain the parameters of our cluster models 
observational data is needed to compare with the models 
predictions. There are a number of observed relations which 
the model should reproduce in order to be an accurate repre- 
sentation of galaxy clusters. For instance the mass tempera- 
ture relation is a fundamental relation which the model will 
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Figure 1 . Core radius plotted against emission weighted temperature for the 8 galaxy clusters and groups 
in the sample of Lloyd-Davics & Ponman (2002) with reliable core radii measurements. The solid line shows 



the relation predicted by the model between core radius and emission weighted temperature when the core 
radius is a constant fraction, 7 percent, of R200 ■ 



need to reproduce. The m ass-temperature relation of Lloyd- _ 
Davies fc Ponman (2002| ) for 20 systems over a large mass 
range was used to compare with the model predictions. The 
masses were derived from analytical models fitted to X-ray 
spectral images and the temperatures are cooling flow cor- 
rected. This is advantageous as our model does not model 
the effects of radiative cooling in cluster cores. 

An observed /3-temperature relation to co mpare against 
our m odel predic tions were obtaine d from Lloyd -DavieE 



Ponman (2002) ( 20 systems) and Helsdon fc P onman 



(2000b| (6 systems). |Lloyd-Davies fc Poiimari (2002| ) fitted 



spherically symmetric gas density profiles to X-ray spectral 
images of galaxy clusters and groups while parameterizing 
either the gas temperature pr ofile or dark matter distribu - 
tion by some functional form. Helsdon fc Ponman (2000b ) 



fitted surface brightness profiles to X-ray images of galaxy 
groups while allowing ellipticity. Their (3 values are th erefore 



only st rictly comparable to the gas density profiles of Lloyd- 



Davies fc Ponman (2002) o n the assumption of isothermal - 
ity. Also the model fitting of Lloyd-Davics fc Ponman (200S ) 
assumes spherical symmetry whereas Helsdon fc Ponman 
(2000b[ do no t, and the temperatures of Lloyd-Davies & 



Ponma 1 (200/ ) are cooling flow corrected whereas those of 



Helsdon & Ponman (2000b) are not. Despite these differ- 



ences the results obtained for these to studies are remarkably 
similar. 

The luminosity-temperature relation data we use to 
compar e wi th ou r model predictions is ta ken from Marke- 
vitch (| 99c]) and Arnaud fc Evrard (1999) for galaxy clus- 



ters and Helsdon fc Ponman (2000a ) for galaxy groups. The 
data of Markevitch (1995) have been corrected for the ef - 
fects of cooling and the sample of Arnaud & Evrard (199£ ) 



by cooling. However the data of Helsdon & Ponman (2000a) 
has not been corrected for the effects of cooling and there 
is not much other data available for these low temperature 
systems. Therefore when comparing them with the model 
predictions we will attempt to quantify the likely effects 
on our results of the galaxy groups being uncorrected. It 
should also be noted that the radii from within which lumi- 
nosity data are extracted are not consistent. This is more of a 
problem for the low mass systems which have flatter surface 
brightness profiles. In these systems the luminosity has gen- 
erally been derived within approximately 0.3 -R200 whereas 
for more massive system the luminosity derived within some 
larger radius. 



4 RESULTS 

A fundamental parameter in all the relations we are inter- 
ested in is the gas temperature. We therefore consider how 
best to extract a gas temperature from our model and then 
go on to compare the predictions of the model with the rela- 
tions listed in Section ^. Finally we will investigate how the 
boundary conditions affect the results in Section 4.5. 



was chosen not to contain any systems significantly affected 



4.1 Temperature derivation 

The effects of raising the specific energy of the gas or increas- 
ing the dark matter concentration on the gas temperature 
are quite similar. In both cases the central temperature and 
temperature gradient are increased. The effects of increas- 
ing the dark matter concentration are more centrally con- 
centrated however, with most of the temperature increase 
occurring near the core of the system. In contrast raising 
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the specific energy of the gas raises the gas temperature out 
to quite large radii. In terms of the effect on the emission- 
weighted temperature increasing the dark matter concentra- 
tion is fairly straight forward, since it has no effect on the 
gas density profile and hence the emission-weighting is not 
significantly altered. Increasing the dark matter concentra- 
tion therefore results in an increase in the emission-weighted 
temperature. In the case of raising the specific energy of 
the gas the flatter gas density profiles result in the emis- 
sion being weighted more towards larger radii where the gas 
temperature is lower. The effect on the emission-weighted 
temperature is therefore not simple and in general depends 
on the radius within which the emission-weighted tempera- 
ture is calculated. For temperatures calculated within -R200 
raising the specific energy of the gas has little or no effect 
whereas temperatures calculated within 0.3 -R200 rise as the 
specific energy of the gas is raised. 

Perhaps the most natural way to derive gas tempera- 
tures from our model is to calculate emission-weighted tem- 
peratures within -R200 since the the observed temperatures 
that we will be using will be weighted in this way. How- 
should be noted that most actual observations of 



ever it 

galaxy bluotoro do not extend to thin radiuo. In the caoo of 
the low mass systems that we are particularly interested in, 
0.3 -R200 is more representative of the radii within which 
emission-weighted temperatures are derived. Even in high 
mass systems the X-ray data rarely extend out to anywhere 
near i?200- We therefore use 0.3 -R200 as the radius within 
which we extract emission-weighted temperatures. It should 
be noted that in general the radius to which X-ray data are 
available increases with increasing system temperature and 
this may lead to some bias in the observed data. As previ- 
ously noted the model does not contain any cooling whereas 
in many observed systems cooling has some effect on their 
emission-weighted temperatures. In some cases it is possible 
to obtain data corrected for the effects of cooling and where 
possible we make use of it. 



4.2 Mass-temperature relation 

An important relation which will constrain our models is 
the mass temperature relation. Raising the specific energy 



of the I pM can affect its temperature, and as the flattening 



of the donoity profilo can puoh gao outoido the virial radiuo, 

the mass within -R200 can also be affected to a certain extent. 
Self-similar theory predicts that the mass should be propor- 
tional toT? and this is ge nerally seen in numerical simula- 
tions (Navarro et al. 1995). Figure |s| shows the total mass 
within -R200 plotted against emission weighted temperature 
within 0.3 -R200 for a range of cluster masses. The systems all 
have a constant concentration parameter of 10. To compare 
our simulated ma ss-temperature relations to obs ervations, 
we use data from Lloyd-Davies & Ponman (2002) for their 



sample of 20 galaxy clusters and groups. The solid line shows 
the mass-temperature relation for the default model. As ex- 
pected from self-similar scaling the relation is a powerlaw 
with an index of 1.5. The other lines show the relations for 
systems where the specific energy has been raised by 0.25 
keV pe r particle (dot-dashed), 0.5 keV per particle (dashed) 



relation is not large, even when the specific energy of the 
gas is raised by a large amount. It should be noted that in 
the 0.75 keV per particle (dotted) case line terminates be- 
fore the lowest mass is reached as beta is flattened to at 
this point. It therefore seems that large deviations from the 
mass-temperature relation due to heating are not possible, 
at least with our model. 

The other process that might affect the mass tempera- 
ture relation is the variation in dark matter concentration. 
The concentration will not affect the mass within -R200 but 
it will have some effect on the temperature. Figure ^| shows 
the total mass within -R200 plotted against emission weighted 
temperature within 0.3 -R200 for several values of the concen- 
tration parameter. It can be seen that as the concentration 
increases from c=2 (solid line) to c=20 (dotted line) the 
normalisation, defined as the mass at 1 keV, decreases from 
1.5 x 10 14 M Q to 3 x 10 13 M Q . This can be compared with 
other relations from observations and theory. The NFW re- 
lation obtained from numerical simulations has a normal- 



isation of 8.7 x 1O"M (Navarro et al. 1995) which is in 



between our two extreme relations. Lloyd-Davies & Pon- 



man (2002) have fitted a T 5 powerlaw to this sample of 



20 galaxy clusters and groups and obtained a normalisation 
of 9.7 ± 1.4 x 10 13 Mq which is also bracketed by our simu- 
lated relations. However this observed normalisation is for all 



galaxy systems. Lloyd-Davies & Ponman (2002) found that 



for systems with temperatures below 4 keV the normalisa- 
tion was 3.1 ± 0.4 x 10 13 Mq which is comparable with our 
high concentration model relations. This can be understood 
in terms of the variation in the concentration parameter with 
system mass. We have seen in Figure |^ that the dark mat- 
ter concentration has a direct effect on the normalisation 
of the mass-temperature relation. Since lower mass systems 
are more concentrated, on average, their mass-temperature 
relation would be expected to have a lower normalisation. 
Whereas rich clusters of galaxies with concentration param- 
eters of around 5 would be expected to have a much higher 
norm alisation to their mass-temperat ure relation, which is 
what Lloyd-Davies & Ponman (2002) observe. The actual 

mean mass-temperature relation should therefore be steeper 

3 

than the simple T 2 relation expected from self-similar scal- 
ing, due to the increase in concentrati on in lower mass sys- 



tems. This effect is indeed observed by Lloyd-Davies & Pon- 



man (2002) who observe the relation to have a logarithmic 



and 0.7 5 keV per particle (dotted). It can be seen that the et al. 200C ; Nevalainen et al 
amount 



slope of 1.96 ± 0.21. 

Figure ^| shows the data with overlayed simulated mass 
temperature relations for different dark matter concentra- 
tion variation prescriptions. The solid line shows the simu- 
lated relation with the concentration parameter varying with 
system mass as specified by Equation |^. It can be seen that 
this relation is much steeper than the self-similar T2 relation 
(dotted line), with a logarithmic slope of ~ 2. The relation 
appears to be a reasonable fit to the observed data, although 
if anything it is slightly steeper. The dashed line shows the 
simulated relation with the conce ntration parameter va rying 
with system mass as predicted by Navarro et al. (1997) for a 
ACDM cosmology (see Equation j). This relation is flatter 
than that predicted by Equation [ , with a logarithmic slope 
of ~ 1.75, s omewhat flatter than is generally observed ( Sato 

200C ; Lloyd-Davies & Ponman 



of deviation from the self-similar mass-temperature 2002). It therefore seems that observations generally favour 
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Figure 2. Total mass within R200 against emission weighted temperature within 0.3 -R200 for simulated 
relations with the specific gas energy raised various amounts; 0.0 keV per particle (solid); 0.25 keV per 
particle (dot-dashed); 0.5 keV per particle (dashed); 0.75 keV per particle (dotted). A constant concen tration 
parameter of 10 is used for all the relations. Observed data (crosses) are taken from the sample of Lloyd 
Davies & Ponman (20021). 
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Figure 3. Total mass within R200 against emission weighted temperature within 0.3 -R200 f° r simulated 
relations with various dark matter concentrations; c=2 (solid); c=5 (dot-dashed); c=10 (dashed); c=20 
( Hotted). The excess energy for a ll the relations is zero. Observed data (crosses) are taken from the sample 
of Lloyd-Davics fc Ponman (2002 ). 
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concentration variation that is steeper than the ACDM pre- 
diction . For this reason the observed relation Equation [| 



therefore derived from the f3 parameters and temperatures 



of the eigh t systems in the sample of Lloyd-Davies & Pon 



will use i to derive the following the results but comparisons 
with the ACDM prediction will be included where necessary. 
The mass-temperature relation appears to be significantly 
affected by the variation in concentration of galaxy clusters 
and these effects should propagate into other relations in- 
volving galaxy cluster temperatures. 

4.3 Beta-temperature relation 

The relationship between the asymptotic slope of gas den- 
sity profile and the temperature of galaxy clusters is im- 
portant since it appears that the main result of raising the 
s pecific energy of the ICM is to flatten the gas density profile 



(Ponman et al. f99£; Lloyd-Davies fc Ponman 2002). Since 



we have previously seen that variation in the dark matter 
concentration of galaxy clusters can affect their emission- 
weighted temperatures, some effect on the beta-temperature 
relation is also to be expected. Figure |E] shows /3 plotted 
against emission weighted temperature within 0.3 -R200 for 
various dark matter concentration parameters. The model 
predictions are shown for concentration parameters of c=10 
(solid), c=15 (dot-dashed), c=20 (dashed) and c=25 (dot- 
ted). In each case the specific energy of the ICM has been 
raised by 0.75 keV per particle. The observed d ata are taken 
from the samples of Helsdon fc Ponman (2000b ) (circles) and 
Lloyd-Davies & Ponman (2002) (crosses). It can be seen that 



as the concentration is increased, the relation is pushed to 
higher temperatures. However the change in the relation is 
not particularly large, and cannot entirely explain the scat- 
ter in ft — T observed. 

The other parameter that can be varied is the amount 
by which the specific energy of the ICM is raised. Figure 
|6| shows (3 plotted against emission weighted temperature 
within 0.3 -R200 with the dark matter concentration vary- 
ing according to Equation ^j. The model predictions for the 

spprifir^ gas pnprgy raiserl by QJ25 kpV per particle (solid) 



0.5 ke y per particle (dot-dashed). 0.75 keV per particle 



(dashed) and 1.0 keV per particle (dotted) are shown. It 
can be seen that there appears to be a considerable amount 
of scatter with none of the energy injection values chosen 
providing a good fit to all the data. In fact the statistical 
errors on the points accounting for only 5 percent of the scat- 
ter about a mean fitted relatio n. This result is not signifi 



cantly changed by removing the Helsdon & Ponman (20001 ) 
points. If this scatter is interpreted as due to a scatter in the 
amount of energy injected into each system then energy in- 
jection covering at the range 0.25 and 1.0 keV per particle 
are required to match the data. Only a small component of 
the scatter can be attributed to scatter in the concentration 
parameters of the systems. 

It is possible to use the predictions of the model to de- 
rive excess energies for individual low mass systems. Since 
the dark halo concentration parameters of the systems affect 
the model predictions to some extent, these measurements 
are best done for system s with known dark matter conc en- 
trations. The sample of Lloyd-Davies fc Ponman (2002j ) is 



the only data available where concentration parameters have 
been measured for low mass systems and also has the ad- 
vantage of the various observational parameters being de- 
rived within a consistent set of radii. Excess energies were 



man (2002 ) with temperatures below 3 keV and are shown in 
Figure [?] plotted against the systems emission weighted tem- 
peratures within 0.3 Jfeoo- It can be seen that apart from the 
three lowest temperature systems the there is little trend in 
excess energy with temperature although there is consider- 
able scatter. The mean excess energy is ~1 keV per particle 
excluding the three lowest temperature points. It is clear 
that the three lowest temperature systems fall considerably 
below this value and show a strong trend of decreasing ex- 
cess energy with temperature. This trend is also noticeable 
in Figures. While it is possible that this is a real effect, the 
role of selection effect in this trend should not be discounted. 

In particular the luminosity of a system decreases as 
excess energy increases and this effect is much more pro- 
nounced in low temperature systems. For this reason if a 
sample is flux limited and there is some distribution in ex- 
cess energies, then there will be a tendency for systems with 
low excess energies to be preferentially selected. The magni- 
tude of this effect will increase with decreasing system tem- 
perature since the range of luminosities for a given range of 
excess energies is larger for lower mass systems. The three 
lowest temperature systems may therefore represent only the 
extreme low end of the distribution in excess energies since 
low temperature systems with higher excess energies would 
not have high enough luminosit ies to be part of the sample of 
Lloyd-Davies & Ponman (2002 ) . It should be noted that one 



of the five higher temperature systems has an excess energy 
higher than the mean by an amount similar to the amount 
that the three low temperature systems are lower than the 
mean. This suggests that there may be enough intrinsic scat- 
ter in the excess energies of the systems to account for this 
effect. 

We can approximate where this flux limit will lie in 
Figure (7] assuming that all these low mass systems are at 
approximately the same redshift and there is some limiting 
luminosity below w hich they cannot be observed. Helsdon 
& Ponman (2000b| ) provide an extensive sample of X-ray 



bright galaxy groups all of which are brighter than 2 x 10 41 
erg s _1 . Using out model a line of constant luminosity in 
the excess energy-temperature plane can be calculated by 
adjusting the excess energy until that luminosity is achieved 
for a range of system masses. This luminosity limit of 2 x 10 41 
erg s" 1 is shown in Figure ^ as the dashed line. It can be 
seen that in general in the sample fall below the line where 
there luminosities should be larger than 2 x 10 41 erg s -1 . 
The drop in excess energy for the lowest mass systems also 
appears to shadow the luminosity limit in the expected way. 
It therefore seems quite likely that this is a selection effect 
rather than a intrinsic property of low mass galaxy groups. 



4.4 Luminosity-temperature relation 

A well studied relation between galaxy cluster properties is 
the luminosity-temperature relation. This is primarily be- 
cause the luminosity and temperature are the two most 
easily measured properties of galaxy clusters. For the high- 
est redshift clusters, which are extremely important from a 
cosmological perspective, these are the only X-ray proper- 
ties that can be measured. It has long been know that the 
luminosity-temperature relation does not follow the predic- 
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Figure 4. Total mass within -R200 against emission weighted temperature within 0.3 -R2011 for simulated 
relations for dark matter concentration varying according to Equation |i| (solid), Equation j| (dashed) and 
a constant c=10 (d otted). The excess energy for al l the relations is zero. Observed data (crosses) are taken 
from the sample of Lloyd-Davies & Ponman (2002). 




Figure 5. Asymptotic slope of the gas density profile, j3, plotted against emission weighted temperature 
within 0.3 R200 f° r simulated relations with various dark matter concentrations; c=10 (solid), c=15 (dot 
dashed) , c=20 (dashe d) and c=25 (dotted). Observed data are taken from the samples of Helsdon & Ponman 



(2000b) (circles) and Lloyd-Davies & Ponman (2002) (crosses). 
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Figure 6. Asymptotic slope of the gas density profile, /3, plotted against emission weighted temperature 
within 0.3 -R200 for simulated relations with specific energies raised to 0.25 keV per particle (solid), 0.5 keV 
per particle (dot-dashed), 0.75 keV per particle (dashed) and 1.0 keV per particle (dotted). In all cases the 
dark matter concentration vari es with mass according to th e relation giv en in Equation Observed data 



are taken from the samples of Helsdon & Ponman (2000b ) (circles) and Lloyd-Davies & Ponman (2005 ) 
(crosses). 
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Figure 7. Excess gas energy within B20Q plotted against emission w eighted temperature within 0.3 -R200 
for the eight systems in the sample of |Lloyd-Davies fc Ponman (2002) with temperatures below 3 keV. The 



dashed line is the model prediction for a constant luminosity of 2 X 10 erg s (see text). 
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tions of self-similar theory. For self-similar clusters, lumi- 
nosity would be expected to be proportional to tempera- 
ture squared (Kaiser 1986). However observational studies 
of the luminosity-temperature relation tend to find steeper 
slopes, with luminosity t ypically proportional t o tempera- 
ture cubed. For instance Edge & Stewart (1991) have mea- 
sure the luminosity-temperature relation for a sample of 



galaxy clusters to be L oc T 2 Ci2±0 10 ; |David et al . (19931) 

r0.297±0.004 



White et al. 



found 



found T oc L v -^' ±u - m and 

T oc L ' 30 * 05 . These results are all in the region of L oc T 3 
It has been suggested that this steepening fr om the predicted 



rela tion might be caused by cool ing flows ( Alien 
1998^1 

c 



Fabian 



Markevitch (1996 ) has derived a luminosity- 



temperature relation corrected for t he effects of cooling flow s 



and measured a slope of 2.63±0.27. Arnaud fc Evrard (199!: ) 
measured the slope of the luminosity-temperature relation 
for a sample galax y clusters without s trong cooling flows to 
be 2.88 ± 0.15 and |Ettori et al. (200l| ) obtain a slope of 2.7 
in the luminosity-temperature relation for a corrected sam- 
ple of galaxy clusters. It therefore seems clear that cooling 
flows cannot provide an explanation for the steepening of 
the luminosity-temperature relation in galaxy clusters. 

Various authors have extended this work down to 
lower mass systems, galaxy groups, and found that the 
luminosity-temperat ure relation steepens in these systems. 



Ponman et al. (1996] ) measured the slope of the luminosity 
temper ature relation for Hickson compact groups to be 



gate the effect of varying the concentration parameter on 
the luminosity-temperature relation, the relations for dark 
matter concentration varying according to Equations ^ and 
^ are shown as the dashed and solid lines respectively. The 
effect of varying the concentration is to significantly steepen 
the luminosity-temperature relation, bringing it into much 
better agreement with the observed data. The relation with 
concentration varying according to Equation ^| has a loga- 
rithmic slope of ~ 2.2 above 2 keV, the one varying accord- 
ing to Equation ^ has a logarithmic slope of ~ 2.7 above 2 
keV. It can be seen that there is some variation in the slope 
with system temperature. The relation with a constant con- 
centration parameter of c=10 (dashed line) has a logarith- 
mic slope of ~ 1.97 although again there is some variation in 
the slope with system temperature. It should be noted that 
strictly the expectation that the luminosity-temperature re- 
lation will have a logarithmic slope of 2 is only true if the 
emission process is pure thermal bremsstrahlung. The con- 
tribution of emission lines at low temperatures will tend to 
flatten the relation so that the relation will be expected only 
to asymptote to a logarithmic slope of 2 at high tempera- 
tures where the effect of emission lines is less important. 

The normalisation of the luminosity-temperature rela- 
tion for varying dark matter concentration in Figure |^ ap- 
pears somewhat higher than the mean trend of the observed 
data. One possible reason for this is that the relat ion we 
use for the core radius in the model is taken from |Lloyd- 



8.2 ± 2.7 while Helsdon fc Ponman (2000b) found the slope Davies fc Ponman (2002j ), whose sample was picked to be 



for a s ample of loose groups to be 4.9 ± O.i 



Helsdon & 



Ponma l (2000a ) have combined a sample of loose and com- 
pact groups to derive a slope of 4.3 ± 0.5. This is sig- 
nificantly steeper than the luminosity-temperature relation 
for rich galaxy clusters. This steepening of the luminosity- 
temperature relation in galaxy groups has generall y been at- 
tribute d to the effect of heating by galaxy winds (Cavaliere 



et al. 1)97; Helsdon fc Ponman 20001 ; Bower et al. 2001) 



There appear to be two separate effects on the luminosity- 
temperature relation which may or may not be related. The 
steeper than predicted slope for high mass systems and the 
steepening of the relation in low mass systems. It is therefore 
extremely interesting to ask how the luminosity-temperature 
relation predicted by our model compares with observations. 
Ho wever in order to calculate the luminosity, some pre- 



as relaxed as possible. It is possible that the value of 7 per- 
cent of R200 which we use is not representative of galaxy 
cluster population as a whole and that less relaxed clusters 
have significantly larger core radii. Previous studies that did 
not limit themselves to relaxed sy stems, have meas ured core 
radii that are considerably larger ( Mohr et al. 1999 ) . Increas- 
ing the core radius for a fixed gas fraction will reduce the 
central gas density and therefore the luminosity. Another re- 
lated point is that the canonical value of j3 = | that we use 
as the default in the model when no heating has occurred 
does not have a great deal of theoretical support. If in fact 
the mean value of (5 in real high mass clusters never quite 
reaches the canonical value and this would also reduce the 
luminosity for a given gas mass. It should be noted t hat the 
mean (5 for the systems a bove 4 keV in the sample of ^loyd- 



scriptioa is needed for the emissivity of a plasma of a given Davies fc Ponman (200S) is 0.61. It is also possible that the 



te mperature, density and m etallicity. We use the method 
of |Knight fc Ponman (1997| ), which uses bi linear interpola- 
tion over the tabulated cooling function of Raymond et al 



(19760to calculate the cluster luminosities and to derive the 
emission weighting for the emission weighted temperatures. 

Figure ^ shows the bolometric luminosity within 0.3 
R200 against emission weighted temperature within 0.3 i?200- 
The dotted line shows the model prediction for systems with 
a fixed dark matter concentration of c=10 and no heating 
of the ICM. This is compared w ith observed data for galaxy 
cluster s from Markevitch (1998) (tria ngles) and Arnaud & 
Evrard | (1999[ ) (squares). The dat a of |Markcvitch (1998|) i s 
corrected for cooling and that of Arnaud $Tl3vrard (199E) 
is selected to contain no systems significantly affected by 
cooling. O bserved data for galaxy gro ups is taken from the 
sample of Helsdon & Ponman (20001) (crosses) and is not 
corrected for cooling. It can be seen that the model relation 
is considerably flatter than the observed data. To investi- 



gas fraction of 0.2 that we are using is somewhat higher than 
the mean value in the systems to which we are comparing 
the model. 

It can be seen in Figure ^ that while the variation of 
the concentration parameter results in the high tempera- 
ture part of the relation being a reasonable fit to the data 
the low temperature part of the relation still over-predicts 
the luminosities of the low temperature systems. To inves- 
tigate whether ICM heating can improve the fit for the 
low temperature systems relations were plotted with vary- 
ing amounts of energy injection into the ICM. Figure ^ 
shows the luminosity-temperature relations for energy in- 
jection of 0.0 keV per particle (solid); 0.25 keV per par- 
ticle (dot-dashed) 0.5 keV per particle (dashed) and 0.75 
keV per particle (dotted). It can be seen that the heating 
the ICM by up to 0.75 keV per particle has little effect on 
the luminosity-temperature relation for the galaxy clusters. 
However the heating does have a significant effect on the 
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Figure 8. Bolometric luminosity within 0.3 -R200 against emission weighted temperature within 0.3 B200 
for simulated relations for dark matter concentration varying according to Equation ^ (solid), Equation 
^ (dashed) and a constant c=10 (dotted). O bserved cooling corrected data are shown for the sample of 
Markevitch (1998) (triangles). The sample of Arnaud & Evrard (1999) (squares) is uncorrect ed but was 



chosen not to c ontain strong cooling flows. The data shown as crosses is for the sample of 
Ponman (2000b) and is not cooling corrected. 



Hclsdon & 



luminosity-temperature relation for the galaxy groups. As 
is to be expected the flattening of the gas density profile 
reduces the luminosity of the low mass systems towards the 
observed relation. It appears that energy injection of be- 
tween 0.25 and 0.5 keV per particle best reproduces the 
observed relation. It should be noted that there is consider- 
able scatter in the observed luminosity-temperature relation 
in galaxy groups, but given the possibility of scatter in both 
the injected energy and concentration parameter this is not 
surprising. 

One discrepancy in this analysis is that the galaxy 
group data is affected to some extent by cooling, while the 
model does not take cooling into account. The luminosity- 
temperature relation is particularly susceptible to the effects 
of cooling since both parameters will be affected by it. We 
have attempted to to quantify the effects of this on our re- 
sults using the observed pro perties of cooling flows in groups. 
Helsdon & Ponman (2002) find that these co ntribute typi - 



cally 2 5 percent of the group luminosity, whilst Lloyd-Davies 
Pom ran (200$ find that removing the effects of cooling 



from spatial-spectral models fitting group X-ray data results 
in an increase in mean temperature of about 20 percent. The 
estimated total average effect of correcting the galaxy group 
data for cooling is shown as an arrow in Figure ^. It can be 
seen that with this correction, most of the data would be 
consistent with a mean injection energy of between 0.5 and 
0.75 keV per particle. 



4.5 Effects of boundary conditions 

The results presented so far are for a specific set of plausible 
boundary conditions, given our present knowledge. The tem- 
perature at -R200 is fixed at 0.57^i r where T V i r is defined by 
Equation and the density at R200 is constant. However it is 
possible that these boundary conditions are over simplified. 
To test the effect of these assumptions the boundary condi- 
tions must be varied in order to see how much difference they 
make to the results. The temperature boundary condition of 
0.5T„; r in particular is based only on numerical simulations 
rather than an empirical evidence. Not only is it possible 
that 0.5 is the wrong normalization but more importantly it 
is possible that it may not scale with T v i r . For example, if 
the intergalactic medium (IGM) has been heated before the 
cluster forms, in the absence of significant cooling the mini- 
mum temperature at -R200 will be the temperature to which 
the IGM was heated. In reality this represents an upper limit 
to the temperature constraint, since if the IGM was heated 
a long time before the system collapsed, the Hubble expan- 
sion would subsequently lower the gas temperature. In low 
mass systems this constraint would result in the tempera- 
ture at 7?2oo having a lower limit of §Ai? where AE in the 
energy injection in keV per particle. To assess how altering 
to the temperature boundary condition might affect the re- 
sults, /3-temperature relations with and without a minimum 
temperature were simulated for an energy injection of 0.75 
keV per particle and dark matter concentration varying ac- 
cording to Equation ^. This is shown in Figure |l(]. The solid 
line shows the relation with no minimum temperature and 
the dashed line shows the effect of imposing a minimum. 
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Figure 9. Bolomctric luminosity within 0.3 -R200 against emission weighted temperature within 0.3 -R200 f° r 
simulated relations with excess gas energy of 0.0 keV per particle (solid); 0.25 keV per particle (dot-dashed) 
0.5 keV per particle (dashed) and 0.75 keV per particle (dotted). In all cases the dark matter concentration 
varies with m ass according to th e relation given in Equation H Observed cooling cor rected data are shown for 
the sample of Markevitch (199£ ) (triangles). The sample of <\rnaud & Evrard (1999 ) (squares) is unc orrected 



but was chosen no t to contain strong cooling flows. The data shown as crosses is for the sample of Helsdon 



Ponman (2000b) and is n ot cooling corrected. The arrow shows is an estimate of the average effect on the 



Helsdon & Ponman (2000k ) data if it was corrected for cooling. 



It can be seen that the result of introducing a minimum 
temperature to the temperature boundary condition is to 
reduce the change in (5 that is needed, to get a specified ex- 
cess energy, for systems below a certain emission-weighted 
temperature. In this case at ~ 2 keV the temperature at 
-R200 hits the limit of |Ai5 and systems below this emission- 
weighted temperature have higher gas temperatures at -R200 
in the raised specific energy case than the default case. This 
causes more energy to go into heating the ICM and less 
into flattening the gas density profile resulting in higher val- 
ues of f3. Increasing the minimum temperature at the outer 
boundary of the model will increase the emission-weighted 
temper ature below which the (3 parameter deviates from the 



predict ons of our standard model. The mechanism by which 



the effect of energy injection on the boundary condition has 
been implemented is very crude, but it does give good indi- 
cation of the general effect that this sort of process would 



have. T he effect is to cause the decrease in (3 with decreasing 



temperature to become less drastic below a certain temper- 
ature. While the number of very low temperature systems 
in our sample is not that large, the data do suggest a steep 
decrease in f3 which flattens off below ~ 1 keV. However 
as noted previously there are likely to be strong selection 
effects at work in the lowest mass system which might be 
expected to cause only the low mass systems with the high- 
est (3 parameters to be observable. 



5 DISCUSSION 

We have modeled two major processes that result in the 
breaking of the expected self-similarity of galaxy clusters; 
energy injection and variation in the dark matter halo con- 
centration. Both these processes are shown to have signifi- 
cant effects on the properties of galaxy clusters. The mass 
temperature relation is not significantly affected by heat- 
ing of the ICM except for the lowest mass systems whose 
temperatures are increased. However variation of the dark 
matter concentration significantly alters the normalisation 
of the mass-temperature relation. T he observed trend in 



concentrat ion parameter with mass (Lloyd-Davies & Pon- 



man 2002), results in a steepening of the mass-temperature 
relation from its expected logarithmic slope of 1.5 to ~ 2. 
This is a good match the observed data, where slopes of this 



order are observed : 1.96 ± 0.21 (Lloyd-Davies fe Ponman 



2002;), 1.79 ±0.1 4 flNevalainen et al. 2000, ) and 2.04 ± 0.04 
( [Sato et al. 2000[ ). The observed mass-temperature relation 



therefore appears to be well matched by the model predic- 
tions when systematic variation in dark matter concentra- 
tion with system mass is included. 

The relation between the asymptotic slope of the gas 
density profile and the emission weighted temperature has 
generally been interpreted as due to energy injection into 
the ICM. The model produces relations that roughly follow 
the trend in the data for energy injection in the range 0.5 to 
TO keV per particle. However the relation has a considerable 
amount of scatter and no one model appears to be consistent 
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Figure 10. Asymptotic slope of the gas density profile, /9, plotted against emission weighted temperature 
within 0.3 -R200 for dark matter concentration varying according to the relation given in Equation ^ and 
raised specific energy of 0.75 keV per particle. The solid line is the relation for the temperature boundary 
condition of 0.5T„i r . The dashed line is the relation when a minimu m temperature of §A_E is im posed on 
the boundary condition. Observed d ata are taken from the samples of Helsdon & Ponman (2000b) (circles) 
and 



Lloyd-Davies & Ponman (2002) (crosses). 



with all the data. This may be due to real scatter in the 
amount of excess energy injected into the systems, although 
it is possible that the errors on the observational data have 
been underestimated. The effect of varying the dark matter 
concentration is to push the relation to higher temperature. 
The inclusion of systematic variation of th e concentration 
parame ter wit h system mass as observed by ( Lloyd-Davies & 
Ponma: 1 2002 ) has the effect of slightly lowering the amount 
of energy that in needed to get a given (5 value. Scatter 
in the concentration of dark matter halos may make some 
contribution to the scatter in the beta-temperature relation. 
However it does not appear to be a large enough effect to 
be the dominant cause. 

The luminosity-temperature relation of the model has 
the expected logarithmic slope of ~ 2 when the model is 
self-similar. However introducing the syst ematic variation 
in dar k matter concentration observed by Lloyd-Davies & 



Ponma: 1 (2002) steepens the slope to 2.7 which is compa- 



rable with the slopes of 2.6 to 2.9 observed by recent au- 



thor s (|Markevitch 199S| ; 



Arnaud & Evrard 1999; Ettori et al 



2001|). This is to be expected as increasing the concentration 
of the low mass systems will increase their emission weighted 
temperatures and so steepen the relation. It therefore seems 
clear that the steepening of the slope of the luminosity- 
temperature relation in the cluster regime is largely, if not 
entirely, explained by the systematic variation of dark mat- 
ter halo concentration with system mass. I t should also be 



noted that the systematic trend seen by Llovd-Davics & X-ray spectr al images of galaxy clusters and groups 



tion of dark matter halos could contribute a considerable 
amount to the scatter in the luminosity-temperature rela- 
tion. However cooling flows are also likely to contribute a 
considerably component to the scatter. 

Models including systematic variation of the dark mat- 
ter concentration with system mass reproduces the slope of 
the luminosity-temperature relation for galaxy clusters but 
does not predict the steepening of the relation seen in galaxy 
groups. However the addition of heating of the ICM steep- 
ens the relation in galaxy groups as flattening the gas den- 
sity profile reduces the luminosity and increases the emission 
weighted temperature. There is considerable scatter in the 
observed relation but most systems appear to fall in between 
the model relations for energy injection of 0.5 keV and 0.75 
keV per particle after correcting for the effects of cooling in 
groups. It should be noted that for galaxy groups, scatter 
in both dark matter concentration and energy injection may 
contribute to the scatter in the relation along with cooling 
flows and other factors. 

It therefore appears that the mean trends in the mass, j3 
and luminosity with system temperature can be explained 
reasonably well by a combination of variation in the dark 
matter halo concentration and energy injection into the 
ICM. The energy injection that best reproduces these trends 
is in the range ~ 0.5 — 0.75 keV per particle. This is com- 
parable with the value o f 0.44 keV per particle found by 
Lloyd-Davies et al. (200C) using analytical models fi tted to 



Ponma: 1 (2002) has considerable scatter which is expected 
theoretically as low mass systems are expected to form a 
higher redshift only on average. Scatter in the concentra- 



Bower 



et al. (2001) derive a value of 0.6-1 keV per particle when 



comparing their model with observational data which is also 
comparable with our result. There is considerable scatter in 
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all the relations which may be due to a combination of scat- 
ter in the concentration parameter and energy injection. 

In the case of the /3-temperature relation it is possible to 
derive excess energies for individual systems and in general 
these are higher than the values derived for the relations as 
a whole. The lowest mass systems however appear to have 
systematically lower excess energies and it is possible that 
this is due to selection effects where low mass systems with 
high excess energies are not observed because they are not 
luminous enough. The mean excess energy of the systems 
below 3 keV excluding the three lowest mass systems is ~1 
keV per particle. This suggests the possibility that care must 
be taken when deriving mean excess energies from cluster 
relations, since the lowest mass systems which tend to have 
the largest influence on the result may have systematically 
lower excess energies due to selection effects. A value of ~ 0.5 
keV per particle therefore appears to be a lower limit on the 
mean excess energy in galaxy clusters and groups, and the 
true mean value may be as high as 1 keV per particle. 

It is possible to estimate the amount of energy that 
could be available as a result of star formation in order to 
assess whether it is possible for this process alone to account 
for the excess energy we observe. Assumi ng that each super- 



nova produces 10 51 erg of kinetic energy ( [Woosley fc Weavei 
1986 here is 0.007 supernova per Mq of stars ( [Bower ct al. 
200 1|) aid the gas and stellar fraction o f galaxy clusters and 



gro ups are 0.2 and 0.11 respectively ( Lloyd-Davies et a^ 
200cj)7}esults in an energy per particle of 1.2 keV. There 
are considerable uncertainties in the assumptions made in 
this derivation but it does suggest that it is possible for star 
formation alone to provide the excess energy we observe, al- 
though the higher mean excess energy of ~1 keV per particle 
we derive for individual systems would require a very high 
efficiency in transferring the energy from supernovae to the 
intergalactic medium. Detailed observation of the metal en- 
richment of the ICM (Finoguenov et al. 2000) also support 
a scenario of heating of the ICM prior to cluster formation 
by galaxy winds. 

It seems clear that the approach of modeling similarity- 
breaking in galaxy clusters and comparing the predictions 
to observations is extremely useful in constraining the pro- 
cesses involved in the formation and evolution of galaxy 
clusters. With the advent of observation of large samples of 
galaxy clusters using XMM-Newton and Chandra it should 
be possible to construct much more detailed models which 
place much greater constraints on the structure of galaxy 
clusters. It may also be possible to directly test the temper- 
ature profiles predicted by our models with such observa- 
tions. 
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